Depinning of a vortex chain in a disordered flow channel 
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We study depinning of vortex chains in channels formed by static, disordered vortex arrays. 
Depinning is governed either by the barrier for defect nucleation or for defect motion, depending on 
whether the chain periodicity is commensurate or incommensurate with the surrounding arrays. We 
analyze the reduction of the gap between these barriers as function of disorder. At large disorder, 
commensurability becomes irrelevant and the pinning force is reduced to a small fraction of the ideal 
shear strength of ordered channels. Implications for experiments on channel devices are discussed. 
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The depinning and dynamics of periodic elastic media 
in a random potential have received a great deal of re- 
cent attention [1]. It was shown, in particular for vortex 
lattices (VL's) in superconductors, that the depinning 
transition in most cases involves plastic deformations [2] . 
A system in which plastic depinning can be studied in 
a controlled way is that of narrow, weak pinning flow 
channels in a superconducting film [3]. In such a system, 
strongly pinned vortices in the channel edges (CE's) pro- 
vide confinement as well as an effective pinning potential 
for chains of vortices inside the channel. By changing the 
magnetic field B, one can vary the ratio between channel 
width and lattice spacing and thus induce incommensu- 
rability between the vortex spacing inside and outside the 
channel. This leads to topological defects in the channel 
which sensitively affect the threshold for plastic flow, as 
evidenced by oscillations of the critical current density 
J c versus field [3]. 

Simulations of channels with perfect hexagonal vortex 
arrays in the CE's, showed sharp peaks in J c for chan- 
nel widths w equal to an integer number n of vortex row 
spacings, i.e. w = nbo- At mismatch (w ^ bo) defects 
occurred in the channel and J c vanished due to the small 
Peierls barrier for defect flow [4] . This however contrasts 
the smooth oscillations found experimentally. Moreover, 
in reality the CE arrays may contain quenched positional 
disorder due to random pinning. This should modify J c 
and its dependence on commensurability. In this Letter 
we consider the most simple near-to-commensurability 
situation w ~ bo, and investigate the effect of disorder 
on depinning of a single vortex chain. We find that com- 
mensurate chains, with a periodicity equal to that of the 
VL in the CE's, depin at a force /„ below the ideal shear 
strength f® by nucleation of defect pairs. At incommen- 
surability, f c is determined by the pinning force of ex- 
isting defects. At weak disorder, a gap is found between 
f n and fd, but for larger disorder commensurability be- 
comes irrelevant and f c saturates at a small fraction of 
/°. This has important consequences for the interpre- 



tation of the experiments [3]. Generally, our results are 
relevant to a wealth of ID problems including models for 
interface growth [5], dynamics of Josephson junctions [6] 
and charge density waves (CDW's). 

We consider a channel at T = with boundaries (CE's) 
formed by two static, disordered arrays with vortex posi- 
tions R„. m = r n>m + d njTn , where r„ !m denote the hexag- 
onal lattice with spacing ao = 2&o/v / 3 (see the inset to 
Fig. 1), the channel width is defined by the spacing be- 
tween rows m = ±1 and d n ^ m are random shifts. We 
restrict ourselves to longitudinal shifts [7] and choose 
d n ,m = d n e x such that the strain (d n +i — d n )/ao is uni- 
formly distributed in the interval [—A, A] with A the 
disorder parameter. The field B sets both the vortex 
density in the CE's (p e = (a 6o) _1 = B/§ with $ the 
flux quantum) and the density p c = (aw) -1 inside the 
channel. Hence, vortices in the channel have an average 
spacing a = aobo/w which can be commensurate with the 
CE arrays (w — bo, a = ao) or incommensurate (a ^ ao). 
The equation of motion for vortex i in the channel is: 



7 a t r, =/-^V7(i 



2 VF(r ( - R„ 



(1) 



V(r) is the interaction potential, j labels other vortices 
in the channel, 7 = B^o/Pf with pf the flux flow resis- 
tivity and / = J$o is the drive along x due to a uniform 
current density J applied perpendicular to the channel. 

For A = the vortex chain can be described by a 
Frenkel-Kontorova model for interacting particles in a 
periodic potential [4]. The ratio between the chain stiff- 
ness and the height of the periodic potential is given by 
g oc A/ao 3> 1 with A the penetration depth. A commen- 
surate chain depins uniformly at a critical force f c = P = 
2aoC66/7i"V / 3 (c66 is the shear modulus) with a velocity 
v = \// 2 — A* 2 /7- At incommensurability, defects of size 
Id = 2Traoy/g S> ao are generated. Their pinning barrier 
and the critical force are essentially vanishing. For f < p 
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and small defect density, given by Cd 
the motion of defects leads to a low mobility regime where 
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v = CdV° d ao with v° d — (-K 2 ^/g/2^)f the pure defect veloc- 
ity [8]. 

We start with a numerical study of the behavior in 
presence of disorder. Eq. (1) was solved using a modi- 
fied London form for V(r) which yields the correct shear 
modulus [4]. We used cyclic boundary conditions, chan- 
nels of length L > lOOOan and we recorded the velocity 
v(f) = (xi)ij and vortex positions Xi(t). 
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FIG. 1. Simulated f-v curves at weak disorder (A = 0.025) 
for w/bo — 1 (•) and w/bo = 0.98 (o) (dashed and full lines 
are the respective results for A = 0). Inset: channel geometry 
with pinned vortices in the gray areas. Their equilibrium po- 
sitions r n ,m are denoted by (o). The disordered arrays (with 
the shifts d x exaggerated for clarity) are denoted by (•). 

The data points in Fig. 1 show f-v curves of a commen- 
surate and an incommensurate chain for A = 0.025. We 
first focus on the commensurate case. Compared to the 
result for ordered CE's (dashed line), f c in presence of 
disorder is clearly reduced. The origin of this reduction 
is that the random strains lower the energy barrier for 
formation of discommensuration pairs in the chain (in- 
terstitial/vacancy pairs in the 2D crystal formed by the 
chain and the CE's). 

We show the depinning process in detail in Fig. 2a 
by plotting the time evolution of the displacements Ui — 
Xi — ia [9]. At t = ti, the force is increased to a value 
/ > ,f c - The motion starts at an unstable site in the 
chain by nuclcation of a vacancy/interstitial pair shown 
as steps of ±an in u. We denote the force at which this 
local nucleation occurs by /„. The process at this site 
repeats periodically with rate R n oc (/ — f n )P and a de- 
pinning exponent fj = 0.46 ±0.04, as previously reported 
for ID periodic media [10]. Due to the nucleation center, 
a domain forms with defect density Cd = R n /{vd) and a 
net velocity v — Cd(vd)ao = -Rn»o with (vd) the average 
defect velocity ~ i>°, see below. In a larger system (Fig. 
2b) , a distribution of unstable sites p with local rate R n:P 
initially leads to the growth of several domains. How- 
ever, when two domains with rates R n ,i > R n ,2 meet, 
their interstitials and vacancies annihilate and domain 1, 



with the larger nucleation rate, then expands at the cost 
of domain 2 with a rate ~ (R\ — R„)- The stationary 
state, covering the entire system, is thus governed by the 
nucleation center with the smallest local threshold /™ m . 
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FIG. 2. Evolution of longitudinal displacements Ui(t) plot- 
ted for clarity in a transverse way vs. x (A = 0.025): (a) 
Nucleation of defect pairs in the commensurate chain of Fig. 
1 (L = 1000a ). (b) Transient response for a system with 
L = 3000ao. The labels 1 and 2 mark competing nucleation 
centers in the long time dynamics. 



Next, we turn to the f-v curves in Fig. 1 at incommen- 
surability (w/bo — 0.98, small defect density Cd ~ l/2Z<j)- 
In contrast to the curve for A = (full line), the data 
for A = 0.025 show a significant threshold force due to 
pinning of the existing defects by random strains in the 
CE's. We define the pinning force for a single defect as 
fd{x) with a distribution {f d } along the channel and a 
maximum value /™ aa: . The shape of the distribution is 
roughly gaussian and has a width ~ 0.1/i. The threshold 
force for a single defect will be f c = f™ ax ■ If the defect 
density is low, the defects occupy only the highest values 
of {f d } and f c < ff ax which in Fig. 1 is f c ~ 0.2/x. 
Once defects are depinned, their velocity Vd fluctuates 
spatially and (vd) is smaller than v®. However, these ef- 
fects strongly decrease with velocity and for / > 2f d nax 
one retains a low mobility regime with dv/df ~ c d v d a , 
as seen in Fig. 1. Considering the regime of larger drive, 
the data exhibits a velocity upturn for / <~ /™". At 
this driving force new defect pairs start to nucleate at a 
strong disorder fluctuation with a rate that exceeds the 
one at which existing defects travel trough the system. 
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FIG. 3. Threshold force f c , defined by a criterion 
v c « O.O2/1/7, versus A for w = f>o (•), w/bo = 0.98 (o) 
and w/bo = 0.95 (A). The thick solid and dashed lines repre- 
sent Eqs. (5) and (6). The typical disorder strengths A c and 
A» are indicated. The insets show f-v curves for A > A c . 

Let us now discuss the disorder dependence of the 
threshold forces. Shown in Fig. 3 arc f c data versus 
A for commensurate and incommensurate chains. For 
w = bg the minimum from many disorder realizations is 
plotted, i.e. f c = f™ m . With growing disorder f™ n 
decreases sharply while for the incommensurate chain 
with w/bo = 0.98, f c <~ f™ ax grows linearly. This be- 
havior changes at A ~ A c where A c is defined as the 
disorder strength where for w = 6 defects first appear 
spontaneously. Above A c both curves (in)decrease more 
slowly with disorder and eventually merge. This change 
in behavior is due to the fact that favorable nucleation 
sites in the commensurate case act as strong pinning site 
for defects in the incommensurate case. The simulations 
show that the presence of pinned defects strongly affects 
the formation of a new nucleus. As a result, /™ m de- 
creases more slowly at larger disorder, i.e. at a higher 
density of pinned defects. The curve for w/bo — 0.95, 
for which the density of static, mismatch induced defects 
is Cd — lj , shows a reduced threshold force and merges 
with the other curves at a disorder strength A = A*. 
For A > A*, disorder induced defects start to overlap 
(i.e. their density becomes > l^ 1 ) and f c decreases fur- 
ther. The large disorder also has a pronounced effect on 
the shape of the f-v curves. As shown in the insets to Fig. 
3, the typical low mobility regime at weak mismatch has 
vanished and all curves exhibit linear behavior, except in 
a small regime for / just above f c . 

To uncover the underlying physics of the phenomena 
described above, we now propose an analytical descrip- 
tion of our system. The energy of a vortex in the channel 
at r = (a;, 0) due to interaction with shifted vortices in 
the CE's is: 

V(r ) = (2tt)- 2 / dkV(k)p e (k)e' kr °. (2) 



V(k) = 27rJ7 /(k 2 + A~ 2 ) with U Q = ^l/2np X 2 , and 
p e (k) are the Fourier transforms of the London poten- 
tial and the vortex density in the edge, respectively. 
For weak disorder (V • d <C 1), p e can be decom- 
posed [11]: p e (r e ,d) ~ (B/$o)(l - V • d + Sp e ) where 
8p e = J2i cos[Ki(r e — d(r e ))] and Kj spans the recipro- 
cal lattice. Substitution in Eq. (2) yields two terms: 
a quasi-periodic potential V p ~ (p/ko) sin[fco(x — d)], 
with ko — 2ir/a [12] due to 5p e of the vortex rows 
nearest to the CE's, and a random, nonlocal contribu- 
tion coming from the density fluctuations: V r (r ) = 
— (B/$ ) / dr e V(r - r e )V • d(r e ) with correlator: 

(V r (0)V r (x)) ~ 4A 2 J7 2 (A/a ) 1+Q e-( a; / A ) 2 . (3) 

Here a depends on the disorder correlations between the 
rows m in the CE's. We choose d x d(x) to be independent 
of row number, which yields a = 2. 
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FIG. 4. (a, b) Mechanical representation of Eq. (4). (a) For 
A < A c a gap exists between the barrier for defect nucleation 
(dashed) and defect pinning, (b) When A > A c disorder in- 
duced defects are always present below f c . The white line 
shows the random phase <j>(x) of the washboard potential. 
(c)top: Evolution of longitudinal displacements of individual 
rows at depinning for w/bo = 3 and A = 0.02. (c)bottom: 
Square lattice representation of the nucleated stacks of dis- 
commensurations. Small arrows indicate the Burgers vector 
of the dislocations terminating each stack. 

To obtain a continuum description of the chain in 
terms of the displacement field u(x), the vortex density 
p c in the channel is decomposed, as was done for p e . 
The resulting interaction of the chain with the CE's is 
h p = Go 1 f( VP + V r )[5 Pc {x,u) - d x u]dx. For A > a , 
both V r and d x u vary slowly while V p and Sp c oscil- 
late rapidly and only two terms in H p remain [13]. The 
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force —SHp/5u on the chain thus contains a commen- 
surate term and a random compression term — d x V r (x), 
which is independent of u. We obtain the equation of 
motion for u by adding the intra-chain force k0 2 u, with 
stiffness k ~ Uon(X/ao) [14], resulting in: 



A, ~ (3V4). 9 - 5/4 , 



(7) 



^d t u = f + nd 2 u — /icos(fc w) — d x V r {x). 



(4) 



The reduced stiffness is now g = «/(fco/ifflo) = 37r(A/ao). 
Eq. (4) describes the transverse displacements of an elas- 
tic string in a tilted washboard potential with random 
phase cj){x) = — dx'V r (x')/n, see Fig. 4. It corre- 
sponds to a commensurate CDW with random forward 
scattering [15] rather than the usual model for a CDW 
in which the commensurability potential is ignored either 
due to strong direct random coupling to u or due to large 
mismatch [16]. 

Next, the dependence of f™ m and f™ ax on the strain 
A can be addressed by considering a deformation in 
the CE's of wavelength Idis'- d x d(r) = A sm(2irx /Idis)- 
When Idis ^> A, the last term in Eq. (4) simplifies to 
-d x V r ~ g 2 AUo/(2V3ldts)cos(27rx/l d i S ). To describe 
nucleation (for / > /u/2), we need the string displace- 
ments S(x) = u — uq relative to the minimum of the tilted 
washboard potential given by u — — [2(1 — f/fj,)] 1 / 2 . 
The restoring force for such displacements is iik 2 uo5(x), 
while the elastic force is nd 2 S. Balancing these with 
d x V r yields S(x) — So cos(2wx /ldi s ). The amplitude of 
the displacement So 
for a wavelength Idis 



Idis/ (I j + k o\ u o\ l dis) is largest 
Id I \/ko\u Q \ resulting in Stf 



V^aog 3 / 2 A/(87ry^o|uo|). Since the critical displacement 
at which the string depins is ~ |«o| [8], the minimum nu- 
cleation force follows from |u | 



xmax. 
°0 



~ 1 - [4A 3 3 / 2 /(5V3) 



4/3 



(5) 



For an existing defect at zero drive, the pinning energy 
is Ed^p(x) = clq 1 j dx'd x Ud(x' — x)V r (x') with Ud the fa- 
miliar shape of a sine-Gordon kink [8] centered at x and 
of size Id- Optimal pinning occurs for deformations in 
the CE's with ldi S — Id, similar to the 'optimal' size of 
a nucleation center for / ~ /x/2. The maximum defect 
pinning force is then given by: 



/r7/« ~ A 3 3 / 2 /V3. 



(6) 



The results (5) and (6) are plotted in Fig. 3, using g = 9 
as in the simulations. Eq. (5) agrees with the numerical 
data, while Eq. (6) can be considered as upper bound. 
The curves merge at A c ~ \/3g~ 3 ' 2 /2, where pinned de- 
fects can appear spontaneously in the channel [17]. 

The disorder strength A* where the density of dis- 
order induced defects becomes ~ l^ 1 can be estimated 
by equating the typical (rms) pin energy of a defect 
with its bare energy ~ (iao^/g. The former is esti- 
mated using the previous form for Ed p and Eq. (3): 
{{Ed^) 2 ) 1 ' 2 ~ 2[/ (3/3^ 2+Q >/ 2 . ( r 1 / 4 A leading to 



also in reasonable agreement with the data in Fig. 3. 
These results show that the effect of disorder rapidly 
grows on increasing the interaction range X/ao- 

Finally, we shortly discuss how these results carry 
over to channels with multiple chains near matching 
(w/bo — n, with n an integer > 2). Without disorder f c 
has sharp peaks at matching of height f® = jibo/w = n/n 
[4]. With disorder, the behavior of f c is similar to that in 
Fig. 3 [18]. In Fig. 4(c) we show the depinning process 
for w/bo = 3, A = 0.02 and / = 0.7/°. As observed, 
the defects involved with depinning consist of stacks of 
discommensurations, coupled between the chains. Each 
stack is terminated by a pair of dislocations with opposite 
Burgers vector along the CE's. This quasi-lD behavior 
at weak disorder can also be described by Eq. 4 by sub- 
stituting ji/n for ji. For larger disorder the threshold 
near matching reduces to ~ 30% of /?. The merging of 
the commensurate and weakly incommensurate thresh- 
old forces now occurs for A™ = A^/n 1 / 4 . For a real- 
istic value of X/ao ~ 4, we estimate A* ~ 0.03 mean- 
ing that rms strains in the CE's < 2% already cause 
strong disorder. At strong disorder and large mismatch 
the quasi-2D nature of the system causes new phenom- 
ena: for w/bo — n± 1/2 transitions from n — > n± 1 rows 
occur, involving pinning of dislocations with misaligned 
Burgers vector [19]. The threshold force then exceeds f c 
near matching [19]. 

In summary, depinning of a vortex chain in a chan- 
nel formed by disordered vortex arrays with (nearly) the 
same periodicity, occurs by nucleation of defect pairs or 
motion of existing defects. The gap between the barriers 
for these two processes and the sharp peak in f c at com- 
mensurability vanish with increasing disorder. For large 
disorder in the channel edges f c saturates at a small frac- 
tion of the ideal lattice strength. 
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